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LONG-TERM  GOAL 

The  long  tenn  goal  of  the  research  project  is  the  understanding  of  the  oceanic  internal  gravity  wave 
field  as  a  balance  of  generation,  transfer  and  dissipation  processes,  and  the  construction  of  a  numerical 
dynamical  model  that  will  predict  changes  of  the  internal  wave  field  in  response  to  changes  in  the 
forcing  and  environmental  fields.  The  internal  wave  model  will  be  based  on  the  integration  of  the 
radiation  balance  equation. 

OBJECTIVES 

The  objectives  of  the  current  project  are  to: 

■  Derive  a  physically  motivated  fonn  of  the  sink  tenn  in  the  radiation  balance  equation  that 
describes  the  dissipation  of  wave  energy  by  wave  breaking,  and 

■  Assess  whether  the  free  parameters  of  this  fonn  can  be  determined  by  reasonably  well-posed 
inverse  problems. 

APPROACH 

The  approach  is  mainly  theoretical  and  numerical.  The  radiation  balance  equation  describes  changes 
of  the  action  density  spectrum  of  the  internal  wave  field  along  wave  group  trajectories  caused  by 
generation,  transfer,  and  dissipation  processes.  The  predicted  quantity  is  the  action  density  spectrum 
as  a  function  of  wavenumber,  position,  and  time.  Because  dissipation  by  wave  breaking  is  a  highly 
nonlinear  process  its  sink  term  in  the  radiation  balance  equation  cannot  be  derived  fully  from  the  basic 
hydrodynamic  equations.  Instead,  such  a  derivation  only  arrives  at  a  functional  form  that  contains  free 
parameters  that  need  to  be  calibrated,  by  comparison  with  observations  or  the  results  of  Large  Eddy 
Simulation  (LES)  models. 

For  the  reflection  off  a  straight  slope  it  can  reasonably  be  assumed  that  dissipation  and  propagation  are 
the  main  processes.  Solutions  of  the  reduced  radiation  balance  equation  can  then  be  constructed  and 
their  parameter  dependence  studied. 
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WORK  COMPLETED 


1 .  Functional  form  of  the  dissipation  sink  term.  Based  on  theoretical  arguments  by  Hassehnann 
(1974)  and  Komen  et  al.  (1994)  a  dissipation  sink  term  of  the  fonn 
Sdiss  =  -  c(k)f  (Ri  -’)  E(k) 

has  been  derived.  Here  k  is  the  wavenumber  vector,  E(k)  the  energy  density  spectrum,  Ri  the 
Richardson  number,  and  c  and / are  two  functions.  The  function  c(k)  detennines  which  wavenumbers 
get  dissipated  by  breaking  events.  We  have  chosen  the  functional  form 

1  if  m  <  m, 


c{k )  =  c 


m  >  m. 


with  three  adjustable  parameters  Co,  mt  and  q,  and  vertical  wavenumber  m.  This  fonnula  reproduces 
the  formula  of  Garrett  and  Gilbert  (1988)  when  q  approaches  infinity.  The  function 
/  (Ri -1)  detennines  the  overall  intensity  of  wave  breaking  and  dissipation.  We  assume/ to  be  a 
monotonic  function  of  Ri  ” 1  and  have  chosen  the  form 


f(m-l)=(Ri  y 


with  an  adjustable  parameter  p.  The  smaller  the  Richardson  number  the  more  vigorous  is  wave 
breaking. 


2.  Solution  of  radiation  balance  equation  for  the  reflection  off  a  straight  slope.  To  determine  the  free 
parameters  c0,  nu,  and  q,p  the  reflection  of  an  incoming  Garrett  and  Munk  internal  wave  spectrum  off 
a  straight  slope  has  been  analyzed.  Because  of  critical  reflection  the  reflected  spectrum  has  (infinitely) 
high  shear  and  zero  Richardson  number  with  associated  vigorous  wave  breaking.  The  dominant 
balance  in  the  radiation  balance  equation  can  therefore  assumed  to  be  between  propagation  and 
dissipation 

vnd  E(k)/dxn  =Sdiss 

where  v  is  the  group  velocity  and  x  the  position.  The  subscript  n  denotes  the  component  nonnal  to  the 
slope.  This  reduced  radiation  balance  equation  has  been  solved  for  different  values  of  the  parameters 
c0 ,  m* ,  q  and  p. 

The  investigation  of  the  four-dimensional  parameter  space  is  facilitated  by  a  certain  separability  of  the 
problem.  First,  one  can  introduce  a  scaled  distance 

i  bg.  z  •  b  g 

which  is  detennined  by  Co  and  p.  The  relation  between  the  energy  E(xn)  (or  any  other  quantity)  and  the 
inverse  Richardson  number  Ri'l(xn)  thus  depends  only  on  m,  and  q.  An  example  is  Figure  1,  which 
shows  E(Ri~l )  for  q= 0  and  q  =  2  and  //  =  m,  b  =50  (b=  1 ,3km). 
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Figure  1:  Energy  density  E  vs.  Ri1  for  q=0  and  q  =2  andu  =  50.  Environmental  parameters  are  f 7 

N=0.06, 

slope  angle  y=0.1,  N  =5.24  * 10 5  s'1.  The  overall  incident  flux  is  17.3  mW  m'2. 

The  dependence  of  Ri  (or  any  other  quantity)  on  xn  is  determined  by  the  parameters  Co  and  p,  where  c0 
determines  the  distance  xlhcrit  where  the  inverse  Richardson  number  decays  to  its  critical  value  and  p 
determines  the  shape  of  this  decay.  Fig.  2  shows  Ri'1  (x„  !  x,hcnt )  for  various  values  of  p. 


Figure  2.  The  inverse  Richardson  number  Ri1  as  a  function  of  x„/X„>critfor  various  values  of  p. 
Environmental  parameters  are  f 7  N  =0.06,  y  =0.1,  N=5.24  * 10 5  s'1.  The  overall  incident  flux 

is  17.3  mWm'2. 


RESULTS 

The  fonn  of  the  dissipation  function  in  the  radiation  balance  equation  needs  to  be  specified  before  any 
meaningful  numerical  experiments  can  be  attempted.  We  have  put  forward  a  theoretically  motivated 
form  of  this  dissipation  function  with  four  free  parameters.  We  considered  the  reflection  problem  off  a 
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straight  slope  where  it  might  reasonable  be  assumed  that  dissipation  is  balanced  by  propagation. 
Analysis  of  the  solutions  as  a  function  of  the  four  parameters  gave  the  following  results: 

1 .  The  total  dissipation  does  not  depend  sensitively  on  the  choice  of  the  parameter  values  and  is 
consistent  with  the  physical  argument  that  all  the  energy  flux  reflected  past  the  “critical  wavenumber” 
is  dissipated.  Figure  3  shows,  that  the  total  dissipation  varies  indeed  by  only  twenty  percent  when  nu 
is  changed  by  more  that  an  order  of  magnitude. 

2.  The  parameters  detennine,  however,  where  in  physical  and  wavenumber  space  the  dissipation 
occurs  and  how  the  resulting  spectra  look  like.  Figure  4  shows  the  vertical  wave  number  spectrum  of 
the  inverse  Richardson  number  at  xn  cnt  for  two  different  values  of  q. 

3.  The  solutions  are  sufficiently  sensitive  to  the  values  of  the  parameters  that  meaningful  inverse 
problems  are  obtained  when  the  theoretical  solutions  are  compared  to  observations. 


Figure  3.  Total  dissipation  Fd  as  a  function  of  m  and  q.  Environmental  parameters  are f 7  N=0.06, 
slope  angle  y=0.1,  N  =5.24  • 10 5  s'1.  The  overall  incident  flux  is  17.3  m  W  m2. 
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Figure  4. 


Vertical  wavenumber  spectrum  of  inverse  Richardson  number  Ri1  at  x,hcrufor  q  =  oo 
(thick  line)  and  q  =  0  (dashed  line).  pp  is  the  vertical  cut-off  wavenumber. 


IMPACT/APPLICATION 

The  development  of  a  predictive  dynamical  model  of  internal  wave  fields  will  have  many  benefits  and 
applications. 

Internal  wave  research  will  benefit  from  such  a  model  since 

■  it  will  provide  understanding  of  the  internal  wave  field  as  a  balance  of  generation,  transfer  and 
dissipation  processes, 

■  it  will  focus  research  (it  is  expected  that  the  proposed  model  will  do  for  internal  wave  dynamics 
what  the  GM  model  did  for  internal  wave  kinematics),  and 

■  it  will  predict  changes  of  the  internal  wave  field  in  response  to  changes  in  the  forcing  and 
environmental  fields. 

The  dynamical  internal  wave  model  can  be  run  in  conjunction  with  circulation  models,  turbulence 
models,  chemical  tracer  models,  and  biological  population  models  where  it  would  predict  the  internal 
wave  induced  transports,  dispersion  and  mixing.  In  conjunction  with  acoustic  transmission  models  the 
model  would  predict  the  internal  wave  induced  “noise.” 
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